SECOND ORDER CORRECTIONS TO MEAN FIELD 
EVOLUTION FOR WEAKLY INTERACTING BOSONS IN THE 
CASE OF 3-BODY INTERACTIONS 
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Abstract. In this paper, we consider the Hamiltonian evolution of N weakly 
interacting Bosons. Assuming triple collisions, its mean field approximation 
is given by a quintic Hartree equation. We construct a second order correc- 
tion to the mean field approximation using a kernel k(t,x,y) and derive an 
evolution equation for k. We show the global existence for the resulting evo- 
lution equation for the correction and establish an apriori estimate comparing 
the approximation to the exact Hamiltonian evolution. Our error estimate 
is global and uniform in time. Comparing with the work of Rodnianski and 
Schlein[22l. and Grillakis, Machedon and Margetis | 12l 1 1 3 1 . where the error 
estimate grows in time, our approximation tracks the exact dynamics for all 
time with an error of the order O (i/Vn\ ■ 



1. Introduction 

In Bosc-Einstcin condensation (BEC), particles of integer spins ("Bosons") oc- 
cupy a macroscopic quantum state often called the "condensate" . The initial ob- 
servation of Einstein and Bose was confirmed experimentally in 1995 and repeated 
later [I] [3l [23] . This phenomenon has stimulated the study of the theory of many- 
body Boson systems. 

In 3d, the dynamics of a system of iV interacting Bosons is governed by a sym- 
metric wave function which solves the N — body Schrodingcr equation 



id t ip N = (A - V N )ip N in 



03^+1 



Or written out explicitly, considering weak interactions and a condensation initial, 
it is 



^ J2 v N {xi-x j ,Xi-x k ) VjvinK 3Ar+1 (1.1) 

i<j<k J 
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in the 3-body interaction case, and 




N 



1 



'^2v N (x i - Xj) ip N in R' 



i<j 



(1.2) 



N 



i=l 



in the 2-body interaction case. The reason that the initial ip N \ t=0 is a tensor 
product for condensate is in p~8j [19l [20] [21] . 

However, as N, the number of particles, increases, solving the N—body Schrodinger 
equation becomes unrealistic. To circumvent this difficulty, we search for effective 
replacements or approximations which are simultaneously easier to approach and 
accurate (in a suitable sense). We are then led to the mean- field approximations 
in which the N — body wave function for the condensate is approximated by ten- 
sor products of a single particle wave function satisfying an appropriate nonlinear 
Schrodinger equation. 

Natural questions arise concerning the justification of the link between the mean- 
field approximations and the actual many-body Hamiltonian evolutions. Elgart, 
Erdos, Schlein, and Yau [4j [5] [6j [7] [10] showed rigorously how mean-field 
limits for Bosons can be extracted in the limit N — > oo by using Bogoliubov-Born- 
Green-Kirkwood-Yvon (BBGKY) hierarchies for reduced density matrices regard- 
ing equation 11.21 By assuming a space-time estimate, Klainerman and Machedon 
gave another proof of the uniqueness part of the argument in |17j . There was sub- 
sequent work by Kirkpatrick, Schlein, and Staffilani [TB] . In Chen and Pavlovic [5] , 
the reduced density matrices approach of equation ll.il was studied. 

Remark 1. It is currently an open question in this field to show directly that the 
limit of the BBGKY hierarchy in 3d satisfies the space-time bound assumed by 
Klainerman and Machedon. The answers in the Id and 2d cases have been found. 
See [2] and [16]. 

Concerning the convergence of the microscopic evolution towards the mean field 
dynamics, Rodnianski and Schlein provided estimates for the rate of convergence in 
the case with Hartrce dynamics by invoking the Fock space formalism of equation 
Oin [22]. 

Recently, following [22], Grillakis, Machedon and Margetis also studied the Fock 
space formalism of equation 11.21 and found a second-order correction (GMM type 
correction) inspired by Wu [24] [25] to the mean- field approximation [12l [13] . This 
paper aims to generalize their result to the Fock space representation of equation 



The main motivation of this paper is to point out that when we apply the GMM 
approximation to the Hamiltonian evolution of many-particle systems equipped 
with 3-body interactions, the error between GMM approximation and the ac- 
tual many-body Hamiltonian evolutions can be controlled uniformly in time. (See 
Knowles and Pickl [T3] for another type of uniform error bound.) We now state our 
problem mathematically. We will discuss the difference bewteen the 2-body and 
3-body cases in subsection 11.31 

1.1. Fock Space Representation of The Problem. First, we set up the Boson 
Fock space T following [12] (T3] [22] . One can also look at Folland [IT] . 
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Definition 1. The Hilbert space Boson Fock space T based on L 2 (R 3 ) contains 
vectors of the form tp = (ip , ip^xi), ip 2 (xi, x%), ■ ■ ■ ) where ip 6 C and ip n £ 
Lg(K 3n ) are symmetric in x±, . . . , x n . The Hilbert space structure of JF is given by 

Definition 2. For f 6 L 2 (R 3 ), we define the (unbounded, closed, densely defined) 
creation operator «*(/) : J~ — > J- and annihilation operator a(f) : J- — >• T by 

1 " 

( a *(/)V>„_l) {xx,X 2 ,--' ,%n) = -7=y2f(Xj)lp n -l(xi,-" ,Xj-l,X j+1 ,"-X n ), 

.7 = 1 



(xi,x 2 ,--- ,x n ) = VnTT J V( n +i)(z,xi, ■ • • ,x n )f(x)dx. 
The operator valued distributions a* and a x are defined by 

a *(f) = / f{x)a* x dx, 



«(/) = / f{x)a x dx. 
These distributions satisfy the canonical commutation relations 

[a>x,ay\ =S(x-y), (1.3) 
[a x ,a v ] = [a*, a*] = 0. 

Define the vacuum state fl € T and the skew-Hcrmitian unbounded operator A 

by 

n = (i,o,o,---) 

A{<j>) = a($)-a*{4>). (1.4) 

It is easy to check that 

e^'O = e-W"/ 8 f 1, ■ • • , (^) V2 • • • 4,(*n), • • • ) • 

This renders the Fock space analogue of the initial data 

N 

v>ati*=o = n^(^) 
i=i 

in eauation ll.il 

In correspondence to equation ll.il we consider the Fock space N particle Hamil- 
tonian Hn : J- —¥ T with the 3-body interaction potential 

H N = J a* x Aa x dx~-^^ J v 3 (x - y,x - z)a* x a* y a* z a x a v a z dxdydz (1.5) 

where the triple interaction is assumed to be symmetric in x, y, and z. As 
shown in [2j, every translation invariant 3-body potential can be written in the 
form vzix — y, x — z). 



Remark 2. Since we use A instead of —A in formula ] 1.51 nonnegative v% corre- 
sponds to the defocusing case. 
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Whence the Fock space representation of equation 1 1.1 1 is the Hamiltonian evolu- 
tion 

e itH Ne -VNA(<f> )^ 

Let the one-particle wave function <fi(t, x) solve the quintic Hartree equation 

d If 22 

i gl^ + A( t > ~ 2^ / v 3K x -y^ x - z )my)\ l<M z )l dydz = (1.6) 

subject to the initial condition 0(0, x) = 4> (x), then the mean field approximation 
for e ltHN e~-' /7rA( -<f>«')Q is the tensor product of <t>{t,x), or 

WMeanField ~ e " 

to be specihe. A derivation of equation [T76] is given in Section [2] 

When the Hamiltonian Hn is the Hamiltonian subject to the two body interac- 
tion 

H N . 2 = J a* x Aa x dx - J v 2 (x - y)a* x a* y a x a y dxdy 

which is the Fock space formalism of equation 11.21 Rodnianski and Schlein [22] 
derived a cubic Hartree equation for <f>(t, x) (equation |1.14[> and showed that the 
mean field approximation works (under suitable assumptions on v) in the sense that 



\Tr 



= 0{ — ) iV^oo; 

where || • \\r r stands for the trace norm in x G R 3 and y G K 3 , and i[> a = e~^ A ^Q. 
For the precise statement of the problem and details of the proof, see Theorem 3.1 of 
Rodnianski and Schlein 22 . Later, in [121 [13], Grillakis, Machedon and Margetis 
introduced a second-order correction (GMM type correction) to the mean field 
approximation of e ltHN ' 2 e~^ A ^°^£l which greatly improved the error. 

Instead of delving into the results in [T^J [13], we state our main theorems first. 
This makes it easier to compare our results with the ones in [T2] [T3J . 

In the main theorem of this paper, we consider the defocusing case: 

v 3 {x -y,x - z) = v(x - y,x - z) 

where 

v(x-y,x- z) = v a (x-y)v (x-z) + v (x-y)v (y- z)+v (x- z)v Q (y-z) (1.7) 

built of a nonnegative regular potential vq which decays fast enough away from the 
origin and has the property that 

v Q (x) = v (-x). 

We remark that our main theorems also work when vq has a proper singularity 
at the origin. To be specific, if for some e G (0, ^) 

«o(*) = frnS > or G 2 +e (z) (1-8) 
\x\ 
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where \ G Co°(R + U{0}) is nonnegative and decreasing and G a the kernel of Bcsscl 
potential, then Theorems [T] holds. Though we currently do not know the physical 
meaning for such potentials if e ^ 0, we would like to understand the analysis 
when singularities appear since the derivation of the quintic NLS uses an interacion 
which goes to a delta function when TV —> oo. Due to the technicality of treating 
the singularities, we restrict to the case of smooth potentials so that the differences 
between the 2-body and 3-body interactions are easier to see. 

Remark 3. For simplicity, let us write A(<j)) as A, A{<j)(t, •)) as A{t), 1)3(0;— y, x—z) 
as v^x— 2,1-3, and 4>{y) as 4> 2 etc. 

1.2. Statement of The Main Theorems. 

Theorem 1. Assume the defocusing case v^,{x — y, x — z) = v(x — y,x — z) where 
v(x — y,x — z) is defined in formula ] 1.7\ 
If <j) , the initial data, satisfies 

(i) finite mass: 

WMli = 1 > 

(ii) finite energy: 

I 2 / 2 2 2 

E ° = 2 / ' V< ^ ' dX+ 6 v(x-y,x- z)\(j) (x)\ \(j) (y)\ |^ ( z )l dxdydz 

sS Ci, 
(Hi) finite variance: 

\\\-\Uli^C 2 , 

then based on the tensor product approximation (mean- field) , we can construct 
ipGMM) a GMM type approximation explained in Theorem^ to the Hamiltonian 
evolution e ltHN e~^ A ^°^Vl for Hjy defined in formula \1.5[ such that the following 
uniform in time error estimate holds 

U GMM - e^e-^Mny <: 



where J- is the Boson Fock space defined in Definition^ and C depends only on v, 
C\ and Ci. 

Remark 4. The construction ofxp GMM does not require v% to have a definite sign. 
However, we take positive sign in Theorem{J\because it leads to a defocusing Hartree 
equation whose global behavior is controllable. 

We prove Theorem [T] via Theorems [5] and [3J stated below. They deal with the 
construction of i^cmm an d the error estimate separately. However, it is worth 
pointing out that Theorem Q] is a special case of Theorems [5] and [31 which apply to 
a more general setting beyond initial data of the form e~^ A ^ ^. 

Theorem 2. Let (f> be a sufficiently smooth solution of the quintic Hartree equation 

i^y> + A0-^Jv 3 (x-y,x-z) |0(y)| 2 \<j>(z)\ 2 dydz = (1.9) 

with initial data <j) Q and the 3-body interaction potential v% being symmetric in x, 
y, and z. Assume the following: 

(1) Let a complex kernel k(t,x,y) € L 2 s (dxdy) for almost allt, solve the equation 

iu t +ug T + gu- (I + p)m = (ip t + [g,p] + um) (I + p)^ 1 u, (1-10) 
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with 

u(t, x, y) := sinh(fc) := k + -^jkkk + ■ ■ • > 
cosh(fc) := I + p(t, x, y) := 5(x — y) + —^kk + . . . , 
g(t, x, y) := -A8(x - y) + (J v 3 (x -y,x-z) \(j)(z)\ 2 dz^j <f>(x)</>(y) 

+ \ (/ v '^ x ~ y,x ~ ^ '^^' 2 \4>{ z )? d v d ^\ s ( x -y)> 

m(t, x, y) := - (J v 3 (x -y,x-z) \(j)(z)\ 2 dz^j <j){x)<j)(y), 

where the products ug T , kk etc. stand for compositions of operators. 

(2) For V defined as in formula \1.5[ the functions 

\\e B Ve- B n\\ Tl \\e B [A, V]e- B n\\ r , \\e B [A, [A, V]]e~ B n\\ r , and \\e B [A, [A, [A, V]]] e - B n\\^ 
are locally integrable in time, where 

B{t) := i J (k(t, x, y)a x a y ~ k(t, x, y)a* x a* y ) dxdy. (1-H) 

(3) J d(t,x,x) dx is also locally integrable in time, where 

d(t, x, y) := (i sinh(fc) t + smh(k)g T + g sinh(fc)) sinh(fc) (1-12) 

— [i cosh(fc) t + [g, cosh(fc)]) cosh(fc) 

— sinh(fc)mcosh(fc) — cosh(fc)msinh(fc). 

Then we define 

$gmm := e -^(*('.-)) e -BW e -i/ t (WxoW+x 1 W)d Sn 

where 

1 

3 



X (t) : = I v 3 (x - y,x - z)\(j)(x)\ 2 |0(y)| 2 \<j){z)\ 2 dxdydz, 



Xi(t) ■ = 7j J d(t,x,x)dx. 
This definition of iPqmm yields the error estimate 

\\1>GMM-e* H "e-^ A ^e- B Wn\\j: 

^ f*\\e B Ve- B nyds f*\\e B lA,V}e- B nyds 
6iV 2 mi 
fi\\e B [A,[A,V]]e- B n\\j:ds fZ\\e B [A,[A,[A,V]]]e- B Sl\\j:ds 
12N 36M 
Theorem 3. Assume v 3 (x — y,x — z) — v(x — y,x — z) i.e. eauation \1.9\ becomes 

iJ^ + A^-i^ j v(x~y,x-z)\(t>{y)\ 2 \<t>(z)\ 2 dydz = {). (1.13) 

If 4> 0l the initial data of quintic Hartree equation ] 1.131 satisfies (i), (ii), and (Hi), 
then the hypotheses in Theorem^ are satisfied globally in time. Moreover, we have 
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the error estimate uniformly in time that 

H G MM-e itH -e-^o) e -Bio )nll:F< C 



where C depends only on v, C\, C% and ||u(0, •, )\\ L 2 

We deduce Theorem [1] from Theorems O and [3] by setting 

k(0,x,y) = 0. 

The proof of Theorem [3] relies on the following theorem regarding the long time 
behavior of the solution to the Hartree equation. 

Theorem 4. If <p solve the Hartree equation \1.13\ subject to (i), (ii), and (in), then 

Ml* <y« fort^l, 
where C is a function of v, C\ and C 2 only. 

1.3. Comparison with Results in [T2J [TB] - In Theorem [21 if we change to 
Hn^2, and equation II .91 to 



ij t <t> + A0 - <j> j v 2 (x - y) \<j>{y)\ 2 dy = 0, (1.14) 

and we let 

g(t, x, y) = -A<5(a; - y) + v 2 {x - y)<j>(t, x)4>(t, y) + (v 2 * |0| 2 )(i, x)S(x - y), 
m{t, x, y) = -v 2 (x - y)4>(t, x)c/>(t, y), 

If 2 2 

Xo(t) = / v z(x - y) \(j)(x)\ \</>{y)\ dxdy, 

then Theorem 1 in reads 

\\ [e -VNA(t) e -B{t) e -if*(N Xo (s)+Xi{s))ds£ l _ ^UH N ,2 e -VNA(0) e -B(0)Q ^ 

^ Jo\\e B V2e- B n\\^d S fi\\e B [A,V 2 ]e- B n\\rds 
N + N i 

If v 2 (x) = X |j r ^' ) , the above error estimate becomes 

|| e -VlV J 4(t) e -S(t) e - l / t (W Xo ( s )+ Xl ( ;s )) ( i Sfi _ e itH N ,2 e -VNA{0) e -B(0) n n ^ ^ 

C(l + t)i 

in [13]. 

Compared with the above long time estimate, Theorem [T] demonstrates that 
there is a substantial difference between the 3-body interaction case and 2-body 
interaction. Technically speaking, the main difference between the 2-body and 3- 
body interactions lies in their error terms. Though the analysis is more involved 
even if we assume smooth potential and the formulas are considerably longer, the 
more complicated error terms in the 3-body interaction case in fact allow more 
room to play. On the one hand, an error term in the 3-body case carries at least a 
pair of u, p or <f> which satisfy some Schrodinger equations, for instance, the term 

\v z (x\ - 7/1,21 - Zi)u(t,x 2 ,xi)u(t, yi,zi)\\ L i L 2 
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in formula 14.21 which can be estimated by Lemma El A typical error term in the 
2-body case can carry only one term of u, p or <f>, for example, the term 

\\vi{xi - yi)u(t, yi,xi)\\ l i ([0 ,t])l2 

implicitly inside formula 47 in [T2]. On the other hand, the error estimate in the 
construction of the second order correction involves L\ and we have no L\ dispersive 
estimates for the Schrodinger equation. Therefore, due to the endpoint Strichartz 
estimates [2], we can construct a Ll(R + ) estimate for the 3-body case which is 
Lemma [5] without having the £2 in the 2-body case which is necessary to apply 
the L\ Kato estimate in [T2l [13] . Or in other words, we do Cauchy-Schwarz in time 
differently. 

For the reason stated above, one can not employ the 3-body case error estimate 
in the 2-body case. Furthermore, the tool of error estimate in the 2-body case 
[121 113) , does not apply to the 3-body case, no matter if V3 is regular or singular 
like formula [1.81 

1.4. Outline of The Paper. We prove Theorem [5] in Section [3] The derivation 
of equations II. 91 and II. 101 is also included there. Section [5] is similar to Sections 3-5 
in [12]. They share the same basic ideas though the computation in this paper is 
more complicated. Therefore we refer the readers to Sections 3-5 in [15] for details 
of the infinitesimal metaplectic representation of symplectic matrices, the main tool 
of Section [21 and the rigorous definition of e B etc. 

Sections [3] and H] are devoted to the proof of Theorem[3] (1) and (3) are verified 
in Section [3] Error estimates are sorted out in Section [4] Like Section 3-4 in |13) . 
we first prove apriori estimates for u which solves equation ll.10) which are needed 
in the error estimates, Section [4] Because the general scheme has been set up in 
[121 [13] , the details of some basic lemmas and theorems are omitted e.g. Theorem 
[5] in Section [3] As mentioned before, our error estimates are done very differently 
from the corresponding ones in [13) . The L\ estimate, Lemma [6] is our key lemma 
for error estimates. 

For smoothness of the presentation of Theorems [5] and [3] we postpone the proof 
of Theorem @] to Section [5] 

2. The Derivation of 2nd Order Corrections / Proof of Theorem [2] 

2.1. Derivation of The Quintic Hartree equation. We first derive the quintic 
Hartree equation II .91 for the one-particle wave function cf> as needed in Theorem [5] 

Lemma 1. The following commutating relations hold, where A denotes A(cj>), and 
A , V are defined by formulas \1.4\ and \1.5\ : 

[A,V] 

= 3 J Va(x — y,x — z)(<f)(x)a*a* z a x a y a z + (j)(x)a* x a* y a* z a y a z )dxdydz 
[A, [A,V]} 

= 6 J v 3 (x - y,x - z)(<j)(x)<j)(y)a*a x a v a z + 2<f>(x)<j>(y)a%a* z a y a z 
+(f){x)(f){y)a x a* y a z a z )dxdydz 
+6 / v^(x — y,x — z) \4>(x)\ 2 a*a* z a y a z dxdydz 
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[A, [A, [A, V}]} 

= 36 J v 3 (x - y,x - z)\(f>(x)\ 2 {<j)(y)a* z a y a z + <j)(y)a*a* z a z )dxdydz 

+6 J v 3 (x - y, x - z){<j){x)<j){y)(f){z)a x a y a z + 4>{x)(j)(y)(t)(z)a* x a*a* z )dxdydz 
+18 J v 3 {x - y, x - z)(cj)(x)(j)(y)(j)(z)a* z a x ay + 4>{x)4>{y)(j){z)a x a y a z )dxdydz 
[A, [A, [A, [A,V]}]} 

= 72 J v 3 (x - y, x - z) \(j){x)\ 2 ((j)(y)(j)(z)a y a z + (j)(y)<j)(z)a* y a z )dxdydz 

+144 J v 3 (x — y,x — z)\<l)(x)\ 2 (j)(y)(j)(z)a* z aydxdydz 

+72 J v 3 (x- y,x- z) \(f>(x)\ 2 \4>(y)\ 2 a* z a z dxdydz 

[A, [A, [A, [A, [A,V]]]]] 
— 360 J v 3 (x — y, x — z) \(/)(x)\ 2 \cf>(y)\ 2 (cf>(z)a z + cf>(z)a* z )dxdydz 

[A, [A, [A, [A, [A, [A,V]]]}]} 
= 720 J v 3 (x - y,x - z)\(j)(x)\ 2 \4>(y)\ 2 \4>{z)\ 2 dxdydz 

Proof. This is a direct calculation using the canonical commutation relation ll.3l □ 

Now, we write * (*) = e^ A{ $e itH *' e-^^e'^Q, for which we carry out 
the calculation in the spirit of equation 3.7 in Rodnianski and Schlein |22) . 

Proposition 1. Let <f> solve the Hartree equation 

d If. 22 

then ^o(t) satisfies 

~*o(t) = (Ho [A, [A, [A, V}}}} - ±N- 2 V - \n-^[A, V] (2.1) 



- ^iV- 1 [A, [A, V}} - ±N~i [A, [A, [A, V}}} 



v 3 (x- y,x-z) \(j>{x)\ 2 \(j)(y)\ 2 \(j>(z)\ 2 dxdydz )* (*) 



N 
+ J 

Proof. Applying the formulas 



(l eC(t) ) ( e ~ c(t) ) = d + h [cA + h i c > [c ' d] 

e c He- c = H+\C,H} + ± [C, [C, H}} + . . . 
to C = yNA and H = H^, we obtain 



~*o(*) =io*o , (2.2) 
i at 
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where 

L Q =-( l e ^ A <A e-^W + e^ A ^H N e-^ A ^ 
i \ot J 

=~ (V/ 2 i + ?L[A, A]) +H + N^IA, H ] + ^[A, [A, H }} - ~ (n~ 2 V + N~ 3 / 2 [A, V] 
+ ^[A, [A,V]] + ^[A, [A, [A, V]}} 
+ I[A, [A, [A, [A, V]}}} + ^{-[A, [A, [A, [A, [A, V]]}]] 
+ ^[A, [A, [A, [A, [A, [A, V}]]]} 



The Hartree equation 11.91 is equivalent to setting terms of order \/N to zero i.e. 
\k + [A, H ] - ~ i [A, [A, [A, [A, [A, V}}}}} = . 
Or more explicitly, the above equation is 

\ j 2 2 I 2 2 

a(i<p t +A(f>--<p / ^3,1-2,1-3 1021 1 03 1 dydz)+a*(i(f> t +A(f>--(f) / u 3) i-2,i-3 |0 2 | 1 03 1 d V dz ) = 

via Lemma[T]and the fact that [A^a^a*] = (AS)(x — y). 
Thus 

t [A, A] + [A, [A, H ]]-~ [A, [A, [A, [A, [A, [A, V]}}}}] = , 
I 5! D 



i.e. equation 12.21 simplifies to 

~*o(t) = (#o - ^[4 [A, [A, [A, V]]]] - iiW- iiV- 3 / 2 [A^] 
- ^N-\A, [A, V}} - ^N~i[A, [A, [A, V}}} 

+ jj « 3 ,l-2,l- 3 |0i| 2 |0 2 | 2 |0 3 | 2 ^2/^)*oW • 

which is equation (2TJ □ 



Because J ^3,1-2,1-3 |0i| 2 1 02 1 2 1 03 1 2 dxdydz only contributes a phase when cj> is 
sufficiently smooth, we write 

^ v?,{x - y,x - z)\(j){x)\ 2 \(j){y)\ 2 \(f){z)\ 2 dxdydz := -N\q 



3 

Then the first two terms on the right-hand side of equation 12.11 are the main ones 
we need to consider, since the next four terms are at most O 

In order to kill the terms involving "only creation operators" i.e. a* a* in 
ji\[A, [A, [A, [A, V}}}}, we introduce B (see ll.lljl and denote 

Hence we have 

i at 

where 
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L = ] T (g e B) e -B + e B Loe -B 



I 



dt 



= L Q - ^N-\ B Ve- B - lN-^e B [A,V]e- B - ^N^e B [A, [A,V]} e - B 

-±N-h B [A, [A, [A,V]]]e- B - N Xo , 
and the quadratic terms 

Lq = \ e- B + e B (h - ±±[A, [A, [A, [A, V]]]]j e~ B . (2.3) 

At this point, we proceed to seek a equation for k s.t. the coefficient of a*a* in 
^ \e B [A, [A, [A, [A, V]}]}e~ B is eliminated. 

Remark 5. One might be concerned of the pure creation a* a* a* in [A, [A, [A, V]]]. 
Lemma\^and the factor will take care of that. Note that [THUS] do not have 

terms like this. 

2.2. Equation for k. 

2.2.1. The infinitesimal metaplectic representation^^. Let sp be the infinite di- 
mensional Lie algebra of matrices of the form 

I d k 
S(d,k,l) = I l _ dT 

where k and I are symmetric, and Quad be the Lie algebra consisting of homoge- 
neous quadratics of the form 



Q(d,k,l) :=~ (a x <) _^ T j ( 



a x a* + a*a x 1 f 

d{x,y) dxdy + - / k(x,y)a x a v dxdy 

- \ J l(x,y)a* x a* y dxdy 

equipped with Poisson bracket. In the spirit of page 185, Folland [11], we define the 
infinitesimal metaplectic representation: a Lie algebra isomorphism I : sp — > Quad 
by Q{d, k, I) = l(S(d, k, I)). Then we see that 

B=T(K) , 



for 



and it follows that 
(i) 

if 1(C) € sp. 
(ii) 



K =Uu° ^ ^'^l, (2.4) 



Jc(t,x,y) 
l(e s Ce- s ) =e I ^l(C)e- x W 



if 1(5) is skcw-Hcrmitian. 



12 



XUWEN CHEN 



(iii) 

e^(a x <) (J) e"^) = (a, a*) e s Q 
if 1(5) is skew-Hermitian. 

Remark 6. Properties (i) and (ii) will be used below, (iii) will be used in Section 


2.2.2. Derivation of Equation \1.1Q[ Use the simplifications noted in Remark [3l 
recall that 

±[A, [A, [A, [A,V]]]} 
= \ / "3,1-2,1-3 |0i| 2 ($2^3 a y a z + <f) 2 (j) 3 a* y a* z )dxdydz 
"3, 1-2, 1-3 |0i| 2 ^^alaydxdydz 



and 



we write 



If 2 2 

2 / U3 ' 1 - 2 > 1 - 3 l^il 1^1 a* z a z dxdydz 
H = J a* x Aa x dx 



n M and M = f - n 
— g J V— m 



with 

.9 = - A(5i_ 2 + (^J "3,1-2,1-3 |0 3 | 2 ^ ~4>l < t } 2 + \ (J "3,1-2,1-3 l</> 2 | 2 1 03 1 2 dydz^ £l-2 

and 



"3,1-2,1-3 I03| 2dz j 0102- 

Of course, we would like to be able to write 

H - I i [A, [A, [A, [A, V]]]] = I (G) + 1 (M) 
Unfortunately, the above equality is not true. For example 

/ -AJ 1 _ 2 \\ y a* Aa x + a^Aa* 
ASr-iJJ-J 2 

However, the commutators of I (G) , 2 (M) and # - ^\[A, [A, [A, [A, V]}]} with B 
are the same as in the discussion in page 287 in [T^j. The same idea applies here. 
Split 

Ho [A, [A, [A,V]]]} = H G +2(M) 

where 

H G = Hq- y « 3 ,i-2,i-3 1^3 1 2 <t>2<t>\ a >l a ydxdydz-^ J "3,1-2,1-3 \4*2? \4> 3 \ 2 a* x a x dxdydz 
which has the property that 

[H G ,B) = [I(G),B], 
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Now, Lq from formula 12.31 reads 

Lq = \^ t e B y- B + e B [H a -^[MMAAA,V]]]]\e- B 
ld_ 



'1 d 



= T \\t. ^e K j e~ K j +H G + (e B H G e- B H G ) + l(e K M e - K ) 

= H G +lffl^- e A e-A + [e B ,H G ]e~ B +X{e K Me~ K ) 

= H G +l(^(ll e K ^ e- K + [e K , G]e~ K + e K Me~ K 
= H G +l{Mi+M 2 +M 3 ) . 

Then by the definition of the isomorphism X, the coefficient of a x a y is — (Mi + 
M 2 +M 3 )i 2 , and the coefficient of a* a* is (Mi+M 2 +M 3 ) 2 i. To write it explicitly: 

- (Mi +M 2 + M 3 )i2 (2.5) 



(Mi +M 2 + M 3 ) 2 i 

(i sinh(fc)t + sinh(fc)g T + g sinh(/c))cosh(/c) 
(zcosh(fc) 4 — [cosh(fc), 5] sinh(fe) 



— sinh(fc)m sinh(fc) — cosh(/c)mcosh(fc) . 
Setting formula [2.51 to confers eauation ll.101 This implies that 

l(Mi+M 2 +M 3 ) = - J(Mi + M 2 + M 3 )u ^ + 2 ^ dxdy 

ixO.* v + a*a x 
d(t, x, y) -— — - — dxdy 

= —J d(t, x,y)a*a x dxdy — — J d(t,x,x)dx 
where d(t,x,y) is given by formula If .121 

Remark 7. (Mi +M 2 +M 3 )ij means the entry on the ith row and the jth column 
of the matrix (Mi + M 2 + M 3 ). 

We summarize the computations we have done so far in this proposition: 

Proposition 2. If (p and k solve equations \1.9\ and \1.10[ then the coefficients of 
a x a y and a x a* in e B [A, [A, [A, [A, V]))]e~ B are and Lq becomes 

Lq =H q - J v 3 (x -y,y-z) \<j)(z)\ 2 <j>(y)(j)(x)a x a y dxdydz 

If 22 
~2 / Vs \ x ~ y->y~ z ) I'M 2 )! a* x a x dxdydz 

1 f 

d(t,x,y)a*a x dxdy — — / d(t,x,x)dx . 
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Recall that * = e B * = e B e VNA ^e itHlf e-^ A (-°h- B ^Sl solves 

±|-* = L*. 

i dt 

We can now write out 

L = L Q -^N- 2 e B Ve- B -^N- 3 / 2 e B [A,V}e- B -^N-'e B [A,[A,V]]e- B 
-L N -k e B [Ai [A) [A,V]]]e~ B - N Xo 

= H Q - J v 3 (x -y,y-z) \(j)(z)\ 2 (j){y)<p{x)a* x a y dxdydz 
1 

~2 



"3(2; -y,y- z) |c/%)| 2 \(j)(z)\ 2 a*a x dxdydz 



-J d(t,x,y)a* y a x dxdy ~ ^N- 2 e B Ve- B - ^N-V 2 e B {A,V] 



e- B 



±N- x e B [A, [A, V]]e' B - ^N~h B [A, [A, [A,V]}}e~ B - ± / ,/(/.,-. ,),/,■ 



= L-Xi~ Nxo 
if we write 

X Q (t) = ~ J v 3 (x- y,x- z)\(f){x)\ 2 \(f){y)\ 2 \(f)(z)\ 2 dxdydz 
Xi(t) = ^ J d(t,x,x)dx. 

Remark 8. Note that (j^j = L and L commutes with functions of time. This is 
needed in the proof of Theorem [H which is below. 

2.2.3. The proof of Theorem^ Applying the above proposition, we can give the 
proof of Theorem [2] at this point. 

^ e -VNA{t) e -B(t) e -i^(Nx (s)+Xi(s))ds n _ ^tH Ne -VNA(0) e -B(0)Q ||^ 
= e i ti(N Xo (s)+Xi(s))ds e B(t) e VNA(t) e tiH N e -VNA(0) e -B(0)Q ||^ 

since e~^ A ^ and e~ Bl ^ are unitary. 
But 

— lift - e l ti {Nxo{s)+Xl{ - s))ds ^ Hi 
dt" 

= 2 Re 



(j /oC^XoM+XiW)*^ _ ^ )e t/a(^Xo(»)+Xi(*))d«^ _ 



= 2 Re (iLft, e < J"o(*xoW+xiW)*'¥ - ft) 
^ 2||£ft ||^|| e </o(*XoM+xi («))*»$ - ft 
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due to the fact that 



Notice that 

m = 

we reach 



-(^ Bv e- B + l N ~ 3/2 ^ A > V ^ B + 
^N~ l e B [A, [A, V]]e- B + ^N~h B [A, [A, [A, V]]} e - B ) o. 



— \m _ e */o( Ar Xo(»)+Xi(s))ds^ II 

dt" 11 
\\e B Ve~ B n\\ r , \\e B [A,V]e- B ny 



W 2 6N? 
|e B [A[A^]]e- S ^lb , \\e B [A,[A,[A,V]}}e- B n\\ r 



12N 

Whence we complete the proof of Theorem[2]because e -^NMt) e -B(t) e -iJ* (N Xo (s)+xx(s))dsQ 
and e^^e-^' 4 ^^-- 8 ' )^ share the same initial data e~^ A ^ e'^tt. 

3. Solving Equation 11.101 / Proof of Theorem [3] (Part I) 

Starting from this section, we begin the proof of Theorem [3] In other words, we 
are assuming that 

v 3 (x - y,x - z) = v(x -y,x- z) 

where v is defined in formula II .71 

We first study equation ll.lOl We prove an apriori estimate for u = sinh(A;) and 
use it in a Duhamel iteration argument to show global existence. Finally we verify 
that J d(t, x, x)dx is locally integrable in time. 

Written in the notations in Remark 02 equation 11.101 reads 

(iu t + ug T + gu — (I + p)m) = (ip t + [g, p] + urn) (I + p) 1 u, 

where 

u(t, x, y) = sinh(fc) = k + T^kkk + • • • , 
cosh(fc)(t, x, y) = I + p(t,x,y) = 81-2 + —^kk + . . . , 

g(t, x, y) = - A5i_ 2 + I / fi-2,1-3 |^> 3 | 2 dz ) 0102 



"1-2,1-3 \(t> 2 \ 2 \03fdydz I Sx-2 



m(t, x,y) = - ij ui_2,i-3 |0 3 | 2 dz J 0102- 

As mentioned in Theorem [2] we write composition of kernels as products in the 
above e.g. 

kk(x 7 y) — / k(x, z)k(z,y)dz. 
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Observe that g(t,x,y) = g(t,y,x), i.e. g* = g; and m(t,x,y) = m(t,y,x), i.e. 



T 

TO = TO. 



Moreover, u T = u, p* = p because k € L 2 s {dxdy). 



Via e K e K = I with X defined in formula 12.41 we obtain the trigonometric 
identity 



uu = cosh(fc)cosh(/c) — / 
= 2p + p 2 
which is a relation between u and p. 

3.1. An Apriori Estimate of u. 

Theorem 5. Let Vs(x — y,x — z) = v(x — y,x — z). If u — sinh(fc) is a solution of 
equation ] 1.1 (A on some time interval [0,T], then there exists a C ^ 0, independent 
ofT, s.t. 

\\u(T)\\ L , <C(1 + K0)|| Z? ). 

The major observation is the following lemma which is also the cornerstone to 
showing Theorem [SJ 

Lemma 2. |13j From equation ! 1 . 1 (A we deduce 

(ip t + [g,p] + urn) (I + pY 1 = -(I + p)' 1 (ip t + [g,p] - mu) (3.1) 
and consequently 

i(uu) t + [g, uu] = mu(I + p) — (I + p)um. (3-2) 

Proof. Multiply equation 11.101 on the right by u, it reads 

(iut + ug T + gu) u— (J + p)niu = (ip t + [g,p] + urn) (I + p)" 1 ^. (3-3) 

Take the adjoint in the operator kernel sense of equation 11.101 multiply on the 
left by u, i.e. 

u {—iu t + g T u + ug) — um(I + p) = uu(I + p)^ 1 {—ipt — [g,p] + TnH) . (3.4) 
Subtracting equations 13.31 and 13.41 we have 

i{uu)t + [g, uu} — (I + p)mu + um(I + p) (3-5) 
= (ipt + [g,p] + um)(I +p)~ 1 uu — uu(I +p) _1 (—ipt — [<?,p] + rri/u) 
With uu = cosh(fc)cosh(fc) — / and uu = 2p + p 2 , we compute 

(I + p)~ l uli - (I + p) = (I + py 1 = uu{I + pY 1 - (I + p) 

and 

(I + pY 1 uu = (I + p)~ 1 p + p = uu(I + p)^ 1 
which transform equation 13.51 to 

i(2p + p 2 ) t + [g, 2p + p 2 } + um(I + p)^ 1 - (I + p) _1 to?Z 
= (ip t + [g, P })((I + p^p + P ) - ((I + py 1 ? + p) (-ip t ~ [g,p]) 

i.e. 

2{ipt + [g,p]) + um(I + pY 1 - (I + p)~ rriu 
= (ip t + [g,p])(I + p)~ 1 p + (I + p)~ 1 p (ip t + [g,p}) 
which is equation 13. II due to I — (I +p)~ 1 p = (I + p) _1 . 
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Multiplying equation 13.11 on the right and left by (I + p) produces 
(I + p)(ip t + [g,p] + urn) = - (ip t + [g,p] - mu) (I + p) 
i.e. equation 13.21 

i{uu) t + [g, uu] = mu{I + p) — {I + p)um 
because uu = 2p + p 2 . 

Taking the trace in formula 13.21 yields 



j^HIz,* = Tr [(l/i)(mu(l+p) - (l+p)ufn)] 



Note that 

Hf 2 = Tr(uu) 



= 2Tr (p) + Tr (p 2 ) 



> M 2 



because p(t, x, y) = njkk + . . . must have a nonnegative trace. So 

^||u||| 2 s$ 2 (Umll^lluH^ + |HIl 2 II u IU 2 IIpIIl 2 ) 

< 2(||m|U a ||u|U a + H| ia ||u||i a ). 
By a Gronwall's inequality, we deduce 

\\u(T)\\ L ^ y) < ^ Umll^^dt+IKO)!!^^ exp (j \\m\\^dt 

The following lemma gives us Theorem [5] 

Lemma 3. If v^{x — y, x — z) = v(x — y, x — z), then 

ll m llL t 1 (R+)L^_ !() ^ C < OO 

Proof. Because 

v(x — y,x — z) 

= v Q (x - y)v (x - z) + v (x - y)v (y - z) + v (x - z)v (y - z), 

we have 

, 2 

2 



HI = J \J v(x-y,x- z)\<t> 3 \ dzj I^J \<t> 2 \ dxdy 

< C [\<j )l \ 2 \cj) 2 \ 2 v 2 (x-y)( fv Q (x-z)\4> 3 \ 2 dz) dxdy 



2 



+C J \<t>i\ 2 \<t> 2 \ v 2 (x-y)U v (y-z)\<p 3 \Uz) dxdy 

+ c J I'/'il 2 !^! 2 (^J - z)v {y ~ z)\(j) 3 \ 2 dz^j dxdy 
I + 11 + III. 
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A combination of Holder and interpolation gives the following estimates 



I + 11 



2C 



< CM* 



v o( x ~ y)\(t>2\ 2 dy 



v l{- - y)\<t>2\ 2 dy 



vq(x — z) |0 3 | dz I dx 



v Q (- - z) |^3 1 dz 



is 



III 



< C>||0||^||0 O ||^||0|| 4 L6 ^C||0 O ||^||0|| 6 L6! 

= C J v (x - zi)v (y - zi)v (x - z 2 )v (y - z 2 ) 
101 1 2 1^2 1 2 I0( z i)| 2 \<t>{z 2 )\ 2 dxdydz x dz 2 

< C dz x dz 2 |0(zi)| 2 \4>{z 2 )\ 2 ( / «o(a: - z\)i^(y - z-Cj^^^dxdy 



Vq( x - z 2)vl{y - z 2 )\<l>- L \ 2 \<l>2\ 2 dxdy 



C\ / dz\4>(z)\* ( / v£(x- z)\4> x \'dx 

2 



fill. 



^0(2; — z^cf)^ dx 



i.e. ||ra||r2 ^ C II 01 1 r © ^ Ci~ 3 , for t ^ 1, by Theorem[4] So we conclude the 
lemma. □ 

Remark 9. Theorem^ also has consequences on p because ||j?||l 2 ^5 ||u||£2. 
3.2. The Existence of u. Because eauation ll.101 

(iu t + ug T + gu — (I + p)m) — (ip t + [g, p] + urn) (I + p) 1 u, 

is fully nonlinear in k, it is not easy to solve for k directly from the equation. 
However, if we put in 

I + p = cosh(fc) = y/I + uu 

in the operator sense, equation 11.101 becomes a quasilinear NLS equation in u — 
sinh(/c). In fact, written out explicitly, the left hand side of equation ll.101 is 



iut + ug + gu 



i-Qi - &x - A,, ) u(t,x,y) 



(3.6) 



v(x-yi,x- z) |0 3 | dz (j>{yi)u(t,y 1 ,y)dy 



+0i 



+02 / u(t,x,yi)[ / v(yi - y,x- z) |<£ 3 | dz)(p(y 1 )dyi 



v.-2,-3 |0 2 | 2 |0 3 | 2 d y dz ) {x)u(t, x, y) 



u— a,— 3 10 2 l 2 103l 2 d y dz ) (yM^ ^ y) 
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and the main term of the right hand side 

d 



ipt + \jg,p] 



at 



p{t,x,y) 



i(x-yi,x- z) |0 3 | 2 dz ) (f>(yi)p(t, yi,y)dy 



(3.7) 



-02 J 2=) yi) fy u (2/i - 2/> a; - z) |0 3 | 2 dz^ 4>(y{)dyi 



W-2,-3|0 2 l I 1 OM^ z, y) 



w— 2,— 3 102 1 I03l dydz) (y)p(t,x,y) 



For our purpose, obtaining some reasonable estimates of u and p — cosh(fc) — I 
is enough. So we would like to get around solving for k and go to u directly. 

But at first, we ask the following: k certainly determines u, but does u determine 
k? The proof of Theorem [2] actually needs a well-defined k. 

We answer the above question by the following lemma: 

Lemma 4. [Tj5] The map 

k i y u = sinh(fc) 

is one to one, onto, continuous, with a continuous inverse, from symmetric Hilbert- 
Schmidt kernels k onto symmetric Hilbert- Schmidt kernels u. 



Proof. The proof of this lemma is in |13) 



□ 



Now we consider the existence of u satisfying equation 11.101 As asserted, equa- 
tion [TTTU] is a quasilinear NLS of u. However, we can transform it into a semilinear 
equation which is easier to deal with, through the following lemma. 

Lemma 5. j!3j The following equations are equivalent for a symmetric, Hilbert- 
Schmidt u: 



iu t + ug + gu = 

T 

iut + ug + gu = 



iu t + ug + gu 



if we set 



W 



F 



(I + p)m + (ip t + [g,p] + um) (I + p) 1 u 
{I + p)ni + — [{I + p) 1 mu + urn (7 + p) 

+ 2 *P* + id,P\, ( i + pV 1 u 
(I + p)rn + — ^(7 + p) 1 mu + urn (7 + p) 1 ^ 

w^i + p)- 1 

/ (uu — z) 1 F (uu — z) 1 VT+ zdz 

2m J 
r 

mu(I + p) — (7 + p)um 



(3.8) 



(3.9) 



77ere, T is a contour enclosing the spectrum of the non-negative Hilbert- Schmidt 
operator uu. 
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Proof. (Sketch) Equation l3.8l is the same as equation ll.101 suitably re-written. The 
keystone of the proof is 

ip t + [g,p] = W. 



But 



ipt + [g,p] = i (Vl + uuj +[g, VT- 



1 

2^i 



(uu — z) (i(uu)t + [g, uu]) (uu — z) \J 1 + zdz 



because 



V / + uu = / (uu — z) V I + zdz 



s-1 



i ((uu — z) ^ + [g, (uu — z) 1 ] = — (uu — z) 1 (i(uu) t + [g, uu]) (uu — z) 

The result follows from equation 13.21 

i(uu)t + [<?, uu] = F = mu(I + p) — (I + p)urn. 

□ 

Whence, we only need to show the existence for equation l3.9l which is of the form 

iu t + ug T + gu = m + N(u) 

where the nonlinear part N(u) involves no derivatives of u. Via the ordinary itera- 
tion procedure, we conclude the following existence theorem: 

Theorem 6. |13j Given uq € L, ^(R 6 ) symmetric, there exists £q such that if 

\\m\\Ll([o,T\)Ll xy) 

then there exists u £ i^°([0, T])L^ x % solving equation \3.9\ and hence equation ] 1 . 1 OA 
with prescribed initial condition u(0, x, y) — uq(x, y) £ L%, ^(R 6 ). 

Since we have shown II "T-ll x, 1 fM.+ li 2 < 00 m Lemma [SJ we can divide R + into 
countably many time intervals [T„,T„ + i] such that ll'Till z, 1 (Tt t , ,])l 2 ^ e o- So 
the above existence theorem in fact implies the global existence of u and thus p. 

Via Theorem O we have 

NU~(R+)£? . < c, 

which implies 

Ibllif (R+)La. iI() < C. 
Moreover, the following estimates hold. 

Theorem 7. Let u E L^°(M. + )L^ x ^ be the solution of equation 1 1 . 1 OA subject to 
uq e L? y )(R ) described in Theorem^ Then u satisfies the following additional 
properties: 

^ - A x - A y ^j u\\ L x t(&+)L 2^ ^ ^ C (3.10) 
i^-A x +A y jp\\ Ll(M+)L 2^ ) < C (3.11) 
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where C only depends on v, C\, C2 and \\uq\\ L 2 ■ See Theorem]]} for C\ and C2. 
Proof. We will only show estimate d. 101 Estimate 13 . 1 II can be shown similarly from 

ipt + [g,p\ = W. 

The proof is separated into 2 parts. 
On the one hand we show 

\\iu t + ug T + gu\\ L i {m+)L ^ ^ ^ ||m|| L i (K+)L ^ ^ + \\N(u)\\ L i {&+)L 2^ ^ C e . 

On the other hand we control the terms in iut+ug T +gu different from (i — A x — A, 
namely 

v{x -yi,x- z) \(f> 3 \ 2 dz ) $(x)(/>{yi)u(yi,y)dyi 



and 



v(x-y,x- z) |0 2 1 2 \cj) 3 \ 2 dydz ) u. 



One sees the above two terms from formula 13.61 
Part I. Recall that 



N(u) = pm + — (jyl + p) 1 mu + um(I + p) ^ u + — W, (I + p) 
We have proven 

\\p\\L°°(W+)L^ y) < ll«lli t -(R+)£?. il0 < Ce- 

Together with the fixed time estimate: 

\\kl\\H-s^\\ k L P \\l\\ H -s (3-12) 

these take care of most of the terms in N(u) because (I + p) 1 and (uu — z) 1 | zS r 
have uniformly bounded operator norms. In inequality 13. 121 ||-||^/_5 stands for the 
Hilbert-Schmidt norm and ||-|| op stands for the operator norm. We only need to 
account for W. However, the fact that \z\ ^ C||u|j? 2 on T implies 



\W\\ Ll(u+)LUy) < C (l+ ||w||| r(R+)i 2 ) \\m\\ Li{M+)Ll y} 



i.e. \\N(u)\\ L i {R+)L 2 <; C. 
Part II. 
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Using Holder, it is not difficult to see the estimate 

J (^J v(x - j/i, x - z) |</> 3 | 2 dz^j 0(x)<p(y 1 )u(y 1 ,y)dy 1 



Ll(R+)Ll 



(x, v ) 





i 

2 ~j 2 


) 4>{x)4>{y 1 )u{yi,y)dyi 





U{x) 1 2 



LJ(R+) 



yi,x- z) |</> 3 | 2 dz I |</>(yi)| 2 dj/i / \u(yi, y)\ 2 dyi I dxdy 



j (yj v(x ~Vi,x- z)\<f> 3 \ 2 dz ^j \<l>( x )\ 2 \0(yi)\ 2 dxd yi I \\u\\ L 2^ 



LJ(R+) 



< a 

It remains to show: 



v) 



J v(x-y,x- z) |0 2 | 2 \<p 3 \ 2 dydz ) a 



C. (3.13) 



Write 



O.y) 



(x - y,x - z) |</> 2 | 2 |0 3 | 2 dj/dz ) dxdy 



y u(x - y,x - z) |0 2 | 2 |0 3 | 2 dydz^j u 
i(t,x,y)\ 2 (^J v(s 
C^J \u{t,x,y)\ 2 (^J v Q (x - y)v (x - z) \(j) 2 \ 2 |0 3 | 2 dydz^j dxdy^j 
+C {^j \u{t,x,y)\ 2 (^j v (x - y)v (y - z) \<p 2 \ 2 \<j) z \ 2 dydz^j dxdy 
+C (J\u(t,x,y)\ 2 (^J v (x - z)v (y - z) \cf) 2 \ 2 \(f) 3 \ 2 dydz^J dxdy 



= I + 11 + III. 
According to the estimate 



J v (x-y)\<t>(y)\^ 



^ c W\\ 2 L6 ^ cr 



we acquire 



I = C ( j \u(t,x,y)\ 2 (^J v a (x-y)\4> 2 \ 2 dy^j dxdy) < Ct i \\u\\ L oo ( j BL+)L 2^^fovt^l, 
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11 + III = 2C I j \u(t,x,y)\ yj v (x - y)v (y - z) \cj> 2 \ \<j> 3 \ dydzj dxdy 

< 2C^J\u(t,x,y)\ 2 ^Jv (x- y )\<b 2 \ 2 dy^ C ||^||* fl dxdyj 

< C^IMlif (R+)i*. ^ for t > 1. 
i.e. estimate 13.131 



v(x -y,x-z) \<j) 2 \ 2 \4> 3 \ 2 dydz ) u 



L\{M + )Ll 



□ 



3.3. The Trace J d(t, x, x) dx. Recall that 

d(t, x,y) = (i sinh(k)t + sinh(fc)g T + g sinh(fc)) sinh(fc) 
— (i cosh(fc) 4 + [g, cosh(fc)]) cosh(fc) 



— sinh(fc)mcosh(fc) — cosh(fc)msinh(fc). 

defined by Formula 1 1.1 21 Rewrite it as 

d(t,x,y) = (iu t + ug T + gu)u-(ip t + [g,p\)(I + p) 
—um(I + p) — (I + p) rriu 

because / commutes with everything and It = 0. 

Notice that if ki(x,y) G L 2 , and k 2 (x,y) G L 2 xy) then 

kik 2 \{x,x)dx = \ / ki(x,y)k 2 (y,x)dy\dx < Pi||l? ||fc 2 || 



At this point, we have already shown that m, iut + ug T + gu, ip t + [g,p] and 
urn G L}(R + )L 2 x . and u, p G L%°(M. + )L 2 X y So except (ip t + [g,p]) I, all traces 
in Formula 1 1.1 21 are well-defined and integrable on IR + . 

However, 



ipt + [g,p] = W, 



for 



W = — ^— [ (uu — z) 1 F (uu — z) 1 vT+ zdz 
2m J 
r 

i* 1 = mu(I + p) — (I + p)ufn. 

Inside the contour integral of W, since (uu— z) 1 | ze r has uniformly bounded 
operator norm and | \/J + 2;| C [l + \\u\\ LOD( ^ + ^ L2 we are in fact dealing 
with 

(Bounded)(H - S){H - S)(Bounded) 

where H—S stands for Hilbert-Schmidt. But (Bounded) (H — S) is Hilbert-Schmidt. 
So we are looking at (H — S)(H — S) which has a trace well-defined and locally 
integrable in time. 
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4. Error Estimates / Proof of Theorem [3] (Part II) 

We finish the proof of Theorem|3]with the proposition below whose proof consists 
of classical techniques. 

Proposition 3. Let 4> to be the solution of the Hartree equation subject to (i), (ii), 
and (Hi). Assume we have 



«^ + A a J 4>\\l,1{m.+)li < C3 , 

d \ 
i-Q t - A* - A y J u\\ L i^+ )L ^ < C* 4 , 



i di~ Ax + Ay ) P ^ L ' (M+)L l,y) ^ ° 5 ' 



d_ 

dt 

then we have the error estimates: 

J \\e B Ve- B fiy dt ^ c 

J \\e B [A,V}e~ B n\\j: dt ^ C 

J \\e B [A, [A,V]]e- B n\\jr dt ^ C 

\\e B [A,[A,[A,V]}}e- B £l\\r dt ^ C 
where C only depends on v, (f>, C3, C4, C5, and 1 1 "Lto | ] ^2 
Remark 10. We can prove 

d_ 

dt 

with the same method to show estimate AS. 1,'ft 

Remark 11. TheoremVA shows that C4, C5 depends only on v, C\, C 2 and ||Mo|| i 2 
So C here is determined by v, C\, C 2 and \\uq\\ L 2 

Remark 12. For Theorem^ we take fc(0,x, y) — i.e. Uq = 0. 

Ideally, we would like to prove Proposition |3] in complete details. However, 



i— + &X I </>|| L l( R +) L 2 < C. 



B * -B B 

e a Xo e = e 



6 Q Tj( ,€. C 



(a x <)(^je- B = (a x <) e K (j) 
{u(xi,xo)a Xa + cosh(fc)(xi,x )a* ^ dxi, 
( a x <) (^j e~ B = (a x a*) e K (^j 
(cosh(k)(x 2 , xo)a X2 + u(x 2 ,x )a* X2 ) dx 2 , 



and 



cosh(k)(x,y) = S(x - y) +p(x,y), 
their products generate a large number of terms. The fact that we will always 
commute the annihilations to the right, e.g. a Xi a V2 a* Z2 = 5(y 2 — z 2 )a* Xi + a x 1 a z 2 a y2i 
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to avoid k(x, x) or related traces, produces even more terms. Hence it is impractical 
to list every single term in e B Ve~ B VL etc., instead, we prove a key lemma and do 
a typical estimate. 

Lemma 6. (Key Lemma) Let xi, 2/1,2/2 G K 3 , £2 G K" 1 ,2/3 G K™ 2 i«if/i i/ie possi- 
bility that n\ or n 2 is zero. Assume f , g satisfy 

d 



Axj. ± A X2 ) f(t,x u x 2 )\\ L i(n+) L 2 sC C, 
d \ 

l gl ± ( A yi+ A V 2 ) ±A y 3 j9(t,yuy2,y3)\\Ll(M+)Ll < C- 



Moreover suppose /|t=o, <?|t=o G L 2 . 
TTien 



di / u (xi - 2/1:^1 - 2/2) \f(t,Xx,x 2 )\ \g(t, 2/1,2/2,2/3)! dxidx 2 dyidy 2 dy 3 ) < C. 



Remark 13. Specializing to the case ni,n 2 = 0,3, or 6, we wz/Z apply Lemma&to 
prove Proposition^ 

In addition to the endpoint Strichartz estimates |14j which are necessary, we 
need the following estimate to prove Lemma [5] 



Claim 1. 



Proof. 



Vq{x - y)vo(x - z)f(y, z)dydz 



L2{S?,dx) 



v o( x - V) v a( x - z )f(y, z)dydz 



vl(x-y) ||w 2 || L 3 11/(2/, OIL* d y 



L?(R 3 ,dx) 



13 



< KLiFoLsii/y =c\\f\ 



We can prove Lemma [6] now. 

Proof. By Duhamel's principle, it suffices to prove 

2 



L~ 



dt 



tt(±A *i ±A ^f(x u x 2 ) 



it(±(A B1 +A V2 )±A 



□ 



ra) 5(2/i,2/2,2/3) 



w (xi-yi,xi~y 2 )dxidx 2 dyidy 2 dy 3 j ^ C \\f\\ L2 \\g\\ L 2 
because we have 

^ ± A ai ± A X2 ± (Aj,! + A V2 ) ± A^J /(i,xi,x2)(7(i,2/i,2/2,2/3)||Li(K+)L 2 x !;) 
< C 

with /|t=o,sit=o G L 2 which also guarantees /,<? G LfL%, by the energy estimate. 
The proof is divided into two steps. 
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Step I: Write the partial Fourier transform to be 
4K) = J e**^f 



then we have 



dx 2 



dx 2 [ d^d^eM^Oe^^ 



« 1 deide a de / a e*' 1 ^-<'O e *(±i)(iex|--|el| >*( fa _£) e *<±WI&l'-|«i| J/fo.fc)/^ 



de 2 



Step II: Let £ 2 > V3 be the phase variables corresponding to x 2 , 2/3- Utilizing Holder 
and Claim [TJ we get 

D 2 (a;i - xi - y2)dxidx 2 dy 1 dy 2 dy 3 )^ 

±u(A yi+ A V2 ) g ^ {yuy2) \l(x x - yi )vl{x l - V2 )dy 1 dy 2 



*S C 



#2 



£ 3 



^ C dt[ d£ 



^ C[ dt d£ 



,±itA« 1 



,±itA„. 



*?3 



e ±lt ( A - +A -)4(2/i,^) 



(yi ^2) , 



dt I dr] 3 



5± «(A w+ A S3 )^ (l/ijW!) 



(yi ,y2 ) , 



^ C II/IIl 2 llfflL 2 (endpoint Strichartz [13]) 

The endpoint Strichartz estimates we used in the last line are the 3d L^L® and the 
6d L^L^. estimates. □ 



4.1. Error term e B Ve B fl, An Example. Write 

'Ve- B 

dx dy dz v(xo - yo,x - z ) 



e B al o e~ B e B a* yo e- B e B a* Zo e~ B e B a Xo e- B e B a Vo e- B e B a Zo e- B (4.1) 



3-BODY SECOND ORDER CORRECTION AND UNIFORM ERROR ESTIMATE 27 

= / ^d^*^^ " 2/0, ,o - zo) 

\u{x\,xo)a Xl + cosh(fc)(xi,xo)a* 1 ^ (u(yi,y )a yi + cosh(k)(yi,yo)a* 1 J 

(u(zi,z )a Zl +cosh(fc)(zi,z )a* 1 ^ (cosh(k)(x 2 , x )a X2 + u(x 2 , x )a* X2 ) 
(cosh(k)(y 2 ,y )a y2 + u{y 2 ,y )a* y2 ) (cosh(fc)(z 2 , z )a Z2 +u(z 2 , z )a* 2 ) 

Because we are applying e B Ve~ B to CI, we neglect the terms in product |4. II which 
have more annihilation operators than creation operators. It is also unnecessary to 
consider terms ending with a Z2 or a X2 a V2 a* Z2 . These facts imply that e B Ve~ B £l has 
nonzero elements solely in its 0th, 2nd, 4th and 6th Fock space slots. To exemplify 
the use of Lemma [51 we estimate two typical terms: the order 6 term 

J dxQdyodzQdxidyidzidx 2 dy 2 dz 2 

v(x - y ,x - z )cosh(k)(xi 1 x )cosh(k)(y 1 ,y Q )cosh(k) 
(z 1 , z )u(x 2 , x )u(y 2 , y )u(z 2 , z )a* Xi a* yi a* Zi a* 2 a* 2 a* 2 

which contributes to the 6th Fock space slot of e B Ve~ B Cl as 
ip 6 {x 1 ,y 1 ,z 1 ,x 2 ,y 2 ,z 2 ) 

= J dx dy dz v(x - y , x - z ) 

cosh(fc)(xi, x )cosh(fc)(yi, y )cosb.(k)(zi, z ) 
u(x 2 ,x )u(y 2 ,y )u(z 2 , z ), 

and an order 4 term 

J dxodyodzodx\dyidzidx 2 dy 2 dz 2 v(xo — yo,xo — zo) 

cosh(fc) (xi , x )cosh(/c) (y t , y )cosh(fc) (zi , z ) u(x 2 , x ) cosh(fc) (y 2 , yo)u(z 2 , z ) 

]|i 4; s|s 4; 

a x 1 a y 1 a z 1 a x 2 a Vi a z 2 

= J dxodyodzodx\dyidzidx 2 dy 2 dz 2 v(xo — yo,xo — zq) 

cosh(fc) (xi , x )cosh(/c) (y! , y )cosh(fc) (z x , z )zt (x 2 , x ) cosh(fc) (y 2 , y )u(z 2 , z ) 
((5(2/2 - Z2)a Xl a yi a Xi a X;i + a Xl a yi a Zi a X2 a Z2 a y2 ) 
which contributes to the 4th Fock space slot of e B Ve~ B Cl as 

ip 4 (xi,yi,z 1 ,x 2 ) = J dx dy dz Q dy 2 v(x - y ,x - z ) 

cosh(fc)(xi, x )cosh(fc)(yi, y )cosh(fc)(zi, z Q ) 
cosh(fc)(y 2 , yo)u(x 2 , x )u(y 2 , z ) 
neglecting symmetrization and normalization. 

4.1.1. Estimate ofip6> a triple product involving one u. Via the fact that 

cosh(fc)(a;, y) = 8(x - y) + p(x, y) 
we write out the product in ip e as 

^6 = ^6,<5<S<5 + V'6,p<5<5 + ^6,pp8 + ^6,ppp 
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according to the factors of cosh carried in each term i.e. 

tpe,888 = J v ( x o -yo,Xo~ Zq)8(xx - x )S(yi - yo)S(zi - z )u(x 2 ,x )u(y 2 ,yo)u(z 2 , z )dxody dz 
= v(xi - yi, xi - zi)u(x2,x 1 )u(y2,yx)u(z 2 , z\) 

and 

ip 6>PP p = J v(x ~yo, x ~z )p(xi, x )p(yi, y )p(zi, z )u(x 2 , x )u(y2, yo)u(z 2 , z )dx a dyodz 
etc. We proceed to estimate the worst term: 

dt (^J \tfj 6 SSS \ 2 dx 1 dy 1 dz 1 dx 2 dy 2 dz 2 

dt ^ \v(xi -y\,xi - z 1 )u(x 2 , xi)u(y 2 , y 1 )u(z 2 , zi)| 2 dx\dy\dz\dx 2 dy 2 dz 2 

< C 

where u(y 2 , yi)u(z 2 , z\) takes the place of g in Lemma El 

For terms in -0 6 involving p, we deal with them as the following: By Cauchy- 
Schwarz on dxodyodzo, we obtain 

\* 

|^6, PP p| dx\dy\dzidx 2 dy 2 dz 2 I dt 
sui>( / \p(xi,x )p(yi,y )p(zi,z )\ 2 dx dyodz dxidyidzi)^ 

i 

\v(x - y ,x - z )u(x 2 ,x )u(y 2 ,y Q )u(z 2 , z )\ 2 dxodyodzodx^y^z^J dt 

where the first integral is majorized by the energy estimate of p, the second integral 
is the same as the one appearing in ip 6 sss and can be taken care of by Lemma [5J 

Remark 14. In the estimate regarding ipQ tPPP , we can do Cauchy-Schwarz in an- 
other way: 

i 

|^6,ppp| dx\dy\dzidx 2 dy 2 dz 2 I dt 
sui>( / \p(x 1 ,x )u(y 2 ,y )u(z 2} z )\ 2 dx dy dz dxidy 2 dz 2 )^ 

\v{x a -yo,x - z )u(x 2 ,x )p(y 1 ,yo)p(z 1 ,z )\ 2 dx dy Q d;:,,<hh d: t d.i ■■, ) dt 

which also works by Lemma\E[ Because \\u\\ ||p|| 
4.1.2. Estimate ofip^, a double product involving one u. 

ip 4 {xi,yi,zi,x 2 ) = dx dy dz dy 2 v(x - yo,x ~ z ) 



cosh(fc) (xi , xo)cosh(fc) (y x , y )cosh(fc) {z x , z a ) 
cosh(fc)(y 2 , yo)u(x 2 , x )u(y 2 , z ) 

4> 4,8888 + ••■ + V'/l.pppp 
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where the worst term is 

ipA.SSSS = J dx dy dz dy 2 v{x ~ yo,xo - z ) (4.2) 

8(xi - x )S(yi - yo)S(zi ~ z )S(y 2 - yo)u(x 2 , x )u(y 2 , z Q ) 
= v(xi - yi,xi - zi)u(x 2 ,xi)u(yi,zi). 
Letting u(yi, z\) be g in Lemma El we derive the desired estimate 

J dt dx 1 dy 1 dz 1 dx 2 \ip 4 j S55 \ 2 ^j < C. 

4.2. Remark for all other error terms. At a glance, we can handle all terms 
using Lemma [HI except 

v(x — y, x — z)<fi(x)(f>(y)(f)(z)a x a y a* z dxdydz 

in [A, [A, [A, V}]}, since all other terms end with a instead of a* . This observation 
allows the application of Lemma [5] But Lemma [5] also applies to 

e B v(x — y, x — z)4>(x)4>(y)cf)(z)a*a*a*dxdydzj e~ B Q. 

because we can let (j>{x\) be f(xx), (j){yi)4>(y 2 ) be g{yi,y 2 ) 

Therefore we have established Proposition [3] and thus Theorem [3] 

5. The Long Time Behavior of The Hartr.ee Equation / Proof of 

Theorem [4] 

In this section, we discuss the Hartrcc equation 11.131 

d If 2 2 

l di^ + 2^ / V<yX ~ V,X ~ Z ' d y dz = ° 

where 

v(x -y,x — z) = v (x - y)v (x - z) + v (x - y)v Q (y - z) + v (x - z)v (y - z). 

We assume the nonnegative regular potential vq decays fast enough away from the 
origin and has the property that 

Vq(x) = Vq(-x). 

Throughout this section, we write 



A = 


J v (x 


- y)vo(x - 


z)\cp{y)\ 2 \4>(z)\ 2 dydz 


B = 


J vo{x 


- y)vo{y - 


z)\<f>(y)\ 2 \0(z)\ 2 dydz 


C = 


J v a (x 


- z)My- 


z)\4>(,y)\ 2 \^z)\ 2 dydz. 



for convenience i.e. equation 11.131 becomes 



ij^ + A0-i(ity-B0-C$=O. (5.1) 



So (ii) becomes 

E(t)\ t=0 = (±J |V0| 2 + l -J{A + B + C) |0| 2 ) | t=0 < oo 
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and (i)-(iii) implies 

E c (t) = ft 2 



^(A + B + C)^ 2 ) <co 



To prove Theorem |4j we are going to argue that 

E c (t) < for t ^ 1 

which leads to 

'E c (t) 



^ c 



V(e l± 



C 



L 2 



t 2 



c 

< — for t ^ 1. 

t 



Here are the details of Theorem |H 



5.1. Conservation of Mass, Momentum, and Energy. First, it is not difficult 
to see the conservation law of the L 2 mass 



where 



and 



dtp - VjP> = 



(5.2) 



because equation 15.11 is of the form 

.8 , 



i—cp + A4> = F{ 



Times —(j) to equation [5T] we acquire 



-po + -(T-Ap+(A + B + C)p = 



where 



Moreover, letting 



Po 



1 

2i 



t - 4>4>t) 



a : = tr(a 3 k) = tr(Vj(f>Vk<t> + Vfc(^Vj( 



A : =(-p + \a+^(A + B + C)p)=Ap-^(A + B + C)p 

e : = l -a+ l -{A + B + C)p 
produces the conservation law of energy 

8 t e — VjCTp + l a = 

where 



(5.3) 



^(A + B + C) Pt -±(A + B + C) t p. 
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A direct computation shows that 



and 



(i)(x') \ 

A t p = / v {x - y)v (x - z)(\cj)(y)\ 2 \<j)(z)\ 2 ) t { — - — )dxdydz 
v (x - y)v (y - z) \(t>{y)\ 2 \(j)(z)\ 2 ( — - — ) t dxdydz 

/I I ^ 

v (x - z)v (y - z) \(j)(y)\ 2 \(t>{z)\ 2 ( — - — ) t dxdydz 



B t p 



v (x - y)v (x - z) \<j>{y)\ \<t>{z)\ ( 



tdxdydz 



v (x - z)v (y - z) \<p(y)\ 2 \(f>(z)\ 2 ( ^ ^ ) t dxdydz 



i.e. 



2 2 10(^)1 

CtP = / «o(a- y)v (y- z) |0(y)| 2 |</>(z)| 2 ( — - — ) t dxdydz 



v (x - y)v (x - z) \4>(y)\ \<f>(z)\ 



2 i „.m2 /I^WI" 



) tdxdydz 



J{A + B + C) t p = 2 J(A + B + C) Pt 
which implies J Iq = i.e. the conservation of energy 

E (t) = l J |V0| 2 + i|(A + S + C)H 2 
Similarly, we derive the conservation law of momentum: 
dtPj - V fe [a] - • /, 

where 

h ■= \{A + B + C) Pj -l(A + B + C) jP . 



(5.4) 



5.2. Conformal Identity. At this point, if we multiply conservation law 15.21 by 

2 



15.41 by ix- 7 and 15.31 by t 2 and add the resulting identities, we obtain the con- 



formal identity: 
where 



d t e c - \7jT j +r = Q 



{ ] -^-)p + tx 3 p j +t 2 e = t 2 



-(A + B + C)p 



V(e I 



r 3 : = {±-±-)p> +tx*a{ + tv> \-t\p+~{A + B + C)p) +t 2 aj, 

2 

r : =t 2 l + txH : j-ntAp + t(n-l)-(A + B + C)p. 



This suggests 



E c + R c = 



(5.5) 
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where 



R c :=t J ({n - + B + C)p + x j l^j dx. 

To determine E c , we calculate 

= I / ^t^ - z )(Pi)jP2P3 - | y x 3 v(x-y,x- z)p 1 (p 2 )jP 3 

= l/ ^3' y ( a; -y l a: -^) a; - 7 [(Pl)i/ J 2-Pl(P2)j] 

= y/ ^^(^ - - ■ v i-2(pip 2 ) 



P 3 «0 - y,a:-«)(a; + y) ■ Vi_ 2 (PiP 2 ) + 3 / PsK^ ~ - y) ■ Vi_ 2 (PiP 2 ) 

-3 / PiP2P3 V i-2w(x- y,a; - 2) • (x + y) - - / p 1 p 2 p 3 V 1 - 2 v(x - y,x - z) ■ (x - y) 

-|« / PlP2P3 v ( x - V: x - z ) 



where Vi_ 2 = V x - y = \{V x - V y ). 
Insert formula [1.71 

v(x —y,x — z) = v (x - y)v {x - z) + v (x - y)v (y - z) + v (x - z)v (y - z) 
to the above computation, it is 

'Ijdx 

"I / Plp2p3V °( X ~ Z ^> ( V i-2 w o(a; - V)) ■ {x + y) - - J p 1 p 2 p 3 v (x - y) (V ' x v (x - z)) ■ (x + y) 
\ J PiP2P3 w o(y - z) (Vi^ 2 v (x ~y))-(x + y) + ^ J p 1 p 2 p 3 v (x - y) (V y v (y - z)) ■ (x + y) 

-| j PiP2P3 v o(y - z) (S7 x v (x - z)) ■ (x + y) + - j p 1 p 2 p 3 v Q (x - z) (V y v Q (y - zj) ■ (x + y) 

2 
3 

-| / PiP2P 3 v o(x - z) (Vi_ 2 uo(a; - y)) ■ (x - y) - - J PiP 2 p 3 v (x - y) (V x v (x - z)) ■ (x - y) 
PiP2P 3 v o{y - z) (Vi- 2 fo(i - y)) ■(x-y) + - J PiP 2 p 3 v a (x - y) (V y v (y - z)) ■ (x - y) 

-g 1 / PiP2P 3 v o{y - z) (\7 x v (x - z)) ■ (x - y) + - J p 1 p 2 p 3 v (x - z) (V a u (y ~ z)) ■ (x - y). 

Notice that, in the above calculation. 

1st + 3rd = 
(4th + 11th) + (6th + 13th) = 

8 f 

(2nd + 9th) + (5th + 12th) = — — / PiP 2 p 3 vo(y — z) (V x vq(x — z j) ■ (x — z). 



n J(A + B + C)p 
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So J xHjdx simplifies to 

xHjdx = -|n J(A + B + C)p 



J- 



-| J PiP2P3 v a( x - z ) ((Vuo) (z - y)) ■ (a; - y) 
\ J PiP2P3 v a{y - z ) ((Vu ) - y)) ■ (a; - ?/) 
■3 / PiP2PsMy - z ) ((Vuo) (a; - 2)) • (a; - z). 



Hence 



i? c = tj(^n- l)^(A + B + C)p + x j l^j dx 

= ~t J p^p^x - z) {v (x -y) + ((Vu ) (a; - y)) • {x - y)} 



~3* / PiP2P3 v o{y - z ) {«o(a; - y) + ((Vt> ) (a; - y)) • (x - y)} 

-|* J PiP2P3 v o(y - z ) {M x - z ) + ((Vv ) (a: - z)) • (a; - z)} . 
When vo decays fast enough, we have 

Rc > 0, 

or in other words 

E c ^ for t ^ 1, 
which implies E c (t) does not increase as claimed. 
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